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Abstract 

The scalar form factor of the nucleon and related physical quantities are in- 
vestigated in the framework of the semibosonized SU(3) Nambu-Jona-Lasinio 
soliton model. We take into account the rotational 1/N C corrections and lin- 
ear m s corrections. The strangeness content of the nucleon in the scalar form 
factor is discussed in detail. In particular, it is found that the m s corrections 
play an essential role of reducing the (iV|ss|iV) arising from the leading or- 
der and rotational 1/N C contributions. We obtain the cr n N(0) = 40.80 MeV, 
Act = a 7TN (2ml) - a nN (0) = 18.18 MeV and (r 2 )^ = 1.50 fm 2 . The re- 
sults are in a remarkable agreement with empirical data analyzed by Gasser, 
Leutwyler, and Sainio Q. 
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I. INTRODUCTION 



Since Cheng 1 l] showed that there is a factor-of-two discrepancy between the empirical 
data for the pion-nucleon sigma term (S^jy) an d the naive estimates of the u-term from 
the mass spectrum, there have been a great deal of discussions and disputes about the H^n 
and a term (see Ref. [fj|| and references therein). Donoghue and Nappi M suggested that 
the discrepancy is due to the presence of strange quarks in the nucleon, i.e. (N\ss\N) ^ 
and showed that (iV|ss|./V) contributes almost 30% to the quark condensate in the nucleon, 
making use of the Skyrme model and bag model. At the first thought, it seems to be 
reasonable, since Cheng used the Zweig rule, i.e. neglected (iV|ss|iV). However, one serious 
question arises: a large fraction of the nucleon mass then stems from strange quarks if one 
follows Ref. 0, which contradicts the quark model. Another assumption was that the ratio 
m s /rh is off by a factor of two, which means that the first order perturbation theory collapses. 
However, this kind of suggestion would lead to a breakdown of the Gell-Mann-Okubo mass 
formula which predict the masses of hadrons in a few percent. 

Motivated by these contradictions, Gasser, Leutwyler and Sainio || recently reanalysed 
the a term prudently, taking advantage of newly accumulated and better nN scattering 
data and considering the strong t-dependence of the scalar form factor a(t) (cr(2m^.) — 
tr(O) ~ 15 MeV). The results of Ref. § were a = 45 ± 8MeV and £ ~ 60MeV. The 
y = 2(N\ss\N) / {N\uu + dd\N), a share of (iV|ss|iV) in the a term, was about 0.2, so 
that the corresponding contribution of the term (jV|ss|iV) to the nucleon mass was about 
130 MeV. 

In the meanwhile, the efforts to understand the a term puzzle theoretically have con- 
tinued [0-0]. However, the bone of contention still lies in the role of strange quarks, more 
specifically the contribution of the (iV|ss|iV) to the a term. Recently, several works insist 
that there is no need to introduce a portion of strange quarks to explain the a term discrep- 
ancy. Bass U proposed that based on the Gribov confinement the value of the a term can be 
explained without need to invoke large strangeness content of the nucleon. Ball, Forte and 
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Tigg || also suggested that with the correct understanding of the baryon matrix element the 
a term (identified with ag = rh(N\uu + dd — 2ss\N) ) can be reproduced without violating 
the Zweig rule. Hence, following these arguments, strange quarks do not contribute to the 
nucleon mass. Though it should be small, it is still important to consider the contribution of 
strange quarks to the a term, in line with recent experiments indicating the fact that strange 
quarks might play an important role of explaining the properties of the nucleon §|||. 

It is the object of the present work to study the strangeness contribution to the a term in 
the framework of the semi-bosonized SU(3) Nambu-Jona-Lasinio soliton model (often called 
as the chiral quark soliton model). In our model, the nucleon is understood explicitly as 
N c valence quarks coupled to the polarized Dirac sea bound by a non-trivial chiral mean 
field configuration. The proper quantum numbers of the nucleon can be acquired by the 
semiclassical quantization |T0|jrT| performed via integrating over the zero-mode fluctuations 



of the pion field around the saddle point. It allows the nucleon to carry proper quantum 
numbers such as spins and isospins. The SU(3) NJL soliton model has a merit in that 
it interpolates between the nonrelativistic naive quark model and the Skyrme model. It 
enables us to study the interplay between these two different models (T2J. The model is 



quite successful in describing the static properties of the baryons and their form factors 

EMI. 

The outline of the paper is as follows: In the next section, we sketch the basic formalism 
for the scalar form factor in SU(3) NJL soliton the model. In section 3, we present the 
numerical results and discuss about them. In section 4, we summarize the present work and 
remark the conclusion. 

II. FORMALISM 

The scalar form factor a(t) is defined as a condensate of u and d quarks in the nucleon: 

a(t) = m(N(p')\uu + dd\N(p)} (1) 
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with m = (jn u + m^)/2 ~ 6 MeV. The t denotes the square of the momentum transfer. Our 
model is characterized by a low-momenta QCD partition function in Euclidean space given 
by the functional integral over pseudoscalar meson and quark fields: 

Z = Jv^V^Vir A exp (- J cPxVUD-Vj (2) 

where 

iD = p{-ip + MW 5 + rh), IT 15 = e inaxa " 15 . (3) 

A a are SU(3) Gell-Mann matrices normalized as TrA a A b = 25 ab . The rh denotes the cur- 
rent quark mass matrix for which we take the form diag(jn u ,m ( i,m s ), where m u , and 
m s are the corresponding current quark masses of the up, down and strange quark, respec- 
tively. Here, we assume that isospin symmetry is not broken, i.e. m u = = rh. The M 
stands for the momentum-dependent dynamical mass arising from the spontaneous chiral 
symmetry breaking. The momentum-dependence of the M introduces the ultra-violet cut- 
off. However, we shall regard it as a constant for simplicity. Instead, we employ a simple 
proper-time regularization. The differential operator iD is expressed in Euclidean space in 
terms of the Euclidean time derivative d T , the Dirac one-particle Hamiltonian H(U) and 
symmetry breaking part [jlT | : 



iD = d T + H(U) + h sb (4) 

with 

a ■ V 

H(U) = —— + (3M U U + /3ml, h sb = /?/i l + /3// 8 A 8 . (5) 

Here, we have made the famous embedding Ansatz for the pseudoscalar fields U 1 ^ and U is 
expressed by 

U = (6) 
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The Uo expresses the SU(2) chiral background field U$ = expi[n ■ rP(r)] with the hedgehog 
Ansatz. P(r) denotes the profile function with proper boundary conditions. fiQ and /i§ are 
defined by /xo = (M a — M u )/3 and /is = — (M a — M u )/\/3. M s and M u are constituent quark 
masses of the s and u quarks respectively. The M u is used as an input parameter, while 
the M s is determined by the gap equation [TTfl. The current strange quark mass m s is also 
settled in the same way. We treat the explicit symmetry breaking term h s b perturbatively. 
The hadronic matrix elements of the 7rN cr-term is related to the correlation function 

a{t) ~ {0 \J N ^l)af N (y,-^)\0) (7) 

at large Euclidean time T. a is the quark operator for the a term, defined by a = m{uu + dd). 



Jn is the nucleon current constructed from iV c quark fields |TIJ 

Jn( X ) — J^ e ii---iN c ^JJ 3 T^i J aiii(x)---' l Pa Nc i Nc (x). (8) 

a,\ • ■ ■ Q!iv c denote spin-flavor indices, while i\ - ■ -in c designate color indices. The matrices 
rjj 3 ^y are taken to endow the corresponding current with the quantum numbers JJ3TT3Y . 
The J B plays the role of creating the baryon state. 

The integral over the quark fields are trivial. The integral over the pseudo-Goldstone 
boson fields can be performed by the saddle point method in the large iV c limit. In order to 
find the quantum 1/N C corrections, it is important to take into account the small oscillations 
of the pseudo-Goldstone bosons around the saddle point and the zero modes. The zero modes 
are taken into account by the soliton expressed by U (x, x 4 ) = A(x 4 )?7(x — Z)A^(x i ) with an 
SU(3) unitary matrix A(t). Hence, the collective action S e ff becomes 

S e// = -iV c Splog(zD) 
= -iV c Splog 



d T + H(U) + A\x A )A{x A ) - idZ ■ V 

(9) 

with the angular velocity A' \xi)A(xi) = iVLe = i£l%\ a /I. Sp denotes the functional trace. 



+ A^x 4 )h sb A{x 4 ) - £(y)/M t (x 4 )-^(v / 2A + A 8 )A(x 4 ) 



The ^ stands for the external scalar field, with regard to which we make a functional deriva- 
tive so as to obtain the sigma form factor: 



a ( t ) = -N c -J--Sphg\d T + H(U) + A\x A )A{x A ) - i{3Z-V 
+ A\x 4 )h sb A{x A ) - £( y )/3A\x 4 )-j=(y/2\ + X 8 )A{xM 



(10) 



It is known that there is the dependence of the a term on the regularization scheme [18 



However, we want to stress the fact that we have employed the proper-time regularization 
and have evaluated possible physical observables such as mass splittings, magnetic moments, 
axial constants and electromagnetic form factors within the same scheme and same values 
of input parameters 1 . Hence, we stick to the proper-time regularization for the a term and 
make use of the same input parameters without adjusting. However, we shall not be here 
bothered by going through all the tedious technical details arising from the regularization 
(see Ref. ]16[| for details). 

Having taken into account the rotational 1/N C corrections and linear m s corrections, we 
arrive at 



a(t)=E su{2) (t)(2 + Dg\A)) N 

+ 4S L /C 1 (t)(Dg ) (A) J R i )^ + ^-K, 2 {t)(D^{A)R p ) 



\/3/i 
4m/i S 

4m/z s 
3^3 



M(t)-£i(t)ir 



Vsi 2 

(D$>{A)DS\A)) 
(D<jL\A)Djg(A)) 



N 



N 



N 



M (t)(D^(A)(D^(A) + l)) N 



im/i 



AAo(t), 



(11) 



where 



Z SU (2)(t) = N c / d 6 x Jo (Qr) 



*l al (x)p* val (x) - Y.\^{En)n{E n )mi{x)^ r 



X 



Ki(t) = yE Jd 3 x Jo (Qr) J d 3 y 



*t ( x )r^ val {x) ■ ^l al {y)i3r^ n {y) 



E n — E, 



val 



x In fact, we have only one free parameter, i.e. the constituent up-quark (down-quark) mass. 
However, it is more or less fixed to around 420 MeV by the mass splitting p3| . 
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+ ~& n (x)T* m (x) ■ ^l{y)(3r^ n {y)n M {E n ,E m ) 



«s(t) = yE Jd 3 xj (Qr) Jd 3 y 



^ m o{x)^ val {x)% al {y)^ m o{y) 



EmO — E, 



vol 



1 



+ ~*i(x)y m o(x)*l (y)p-* n (y)ll M (E n ,E m 



• A/ 'i(t) = ^E / d 3 xj (Qr) d 3 y 



mi(x)(3rm val {x) ■ ^l al (y)f3T^ n (y) 



E n E va i 



+ -#t(ar)/3f# m (z) . ^l n {y)(3r^ n {y)np{E n ,E n 



^(t) = ^J2 Jd 3 xUQr) Jd 3 y 

n,m° 



^Lo {X)^val (X) Wgm° (y) 

E m o — E val 



-^Ux)f3^ m o(x)¥ m0 (y)^ n (y)n p (E n ,E m 



Ao(t) 



3iV, 



E va i 



1 



+ -^(x)/^ m (x) ■ *l(y)P* n (y)Tlf,(E n ,E„ 



(12) 



The subscripts i and p in the collective part are i = 1,2,3 and p — 4,5, 6, 7, respectively. 
Ii and Ki are respectively the moments of inertia and anomalous moments of inertia ||13|| . 
When t — > 0, /Cj(£) become K^. The £5(7(2) corresponds to the 7iN sigma term in SU(2) [I3| 
at £ = 0, which can be obtained by the Feynman-Hellman theorem 

dE(fh) 



^SU(2) 



m- 



dm 



(13) 



m=0 



where E stands for the classical soliton energy The regularization functions R(E n ), 
R M (E n , E m ), Rp(E n , E m ) are defined by 

du 



K(E n ) 
T^MyEn: E m ) 



Tlf3{E n , E n 



(u;ki)\E n \e 



-uEi 



1 sign(£ n ) - sign(£ m ) 



E n E m 



00 du 



2\fKU 



4>{u;Ai) 



pp p -uEl _ rp -uEl 

E n E m 



(14) 



2 R_M(E n , E m ) is not actually a regularization function, since /Q come from the imaginary part of 
the action. It does not depend on the cut-off parameter. 
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respectively. The () n stands for the expectation value of the Wigner D functions in collective 
space apanned by A. The expectation values of the D functions can be evaluated by SU(3) 
Clebsch-Gordan coefficients found in Refs. [p]| , pO[| . With SU(3) symmetry explicitly broken 
by m s , the collective part is no longer SU(3)-symmetric. Therefore, the eigenstates of the 
hamiltonian are not in a pure octet or decuplet but mixed states. Since we treat the strange 
quark mass m s perturbatively, we can obtain the mixed SU(3) baryonic states as follows: 

|8,iV> = |8,iV) + 4|10,iV) + 4\27,N) (15) 

with 

cf Q = ^(o- - ri)/ 2 m s ,c^. = ^(3a + r 1 - 4r 2 )I 2 m s . (16) 

The constant a is related to the £5(7(2) by T,su(2) = 2/3(m u + m^a. r« denotes the ratio 
Ki/U 

Since the Cheng-Dashen point is out of the physical region, it is necessary to extrapolate 
to the region t > 0. This can be done by the analytic continuation of the \q\, i.e. \q\ — > 
i\q\ so that we may have the positive t up to the Cheng-Dashen point (t = 2m\). The 
analytic continuation above the threshold t = is not valid in our model, since above 
this threshold, the correlation between mesonic clouds is getting important p2 . Hence, 



in this work, we only evaluate the scalar form factor from the Cheng-Dashen point to the 
physical channel (space-like region: t < 0). 

III. NUMERICAL RESULTS AND DICUSSION 

In order to calculate the <J n N{t) numerically, we take advantage of the Kahana-Ripka 



discretized basis [ 2fJ . Figure 1 shows the scalar form factor as a function of the constituent 
quark mass M = M u = M^. The a(t) decreases as the M increases, in particular, below 
t = 0. As a result, the difference between the a(2m^) and cr(0) changes drastically when we 
increase the M from 370 MeV to 450 MeV, as shown in Table 1. We select the M = 420 MeV 



S 



for the best fit as we did for other observables. The error bar presented in Fig. 1 stands 
for the empirical analysis due to Gasser, Leutwyler, and Sainio ||, i.e. cr(0) = 45 ± 8 MeV. 
Our numerical prediction is in a remarkable agreement with Ref. ||. It is also interesting 
to see how the m s corrections contribute to the scalar form factor. As shown in Fig. 2, 
the m s corrections are very small. At t — 0, the m s corrections contribute to the a term 
about 2% which is negligible. However, the m s corrections play a significant role of reducing 
remarkably the large strangeness contribution (JV|ss|iV) arising from the leading term and 
rotational 1/N C corrections. With the m s corrections taken into account, we obtain y = 0.27 
in case of the M = 420MeV, which agrees with the empirical value y ~ 0.2 || within about 
30%, whereas we have y = 0.48 without the m s corrections. It is already known that the 
explicit symmetry breaking term quenchs the (N\ss\N) [p?) - 29 1 . 

The difference Aa = <j(2m 2 T ) — <r(0) we have obtained is 18.18MeV. This value is very 



close to what Gasser and Leutwyler extracted ||21|| , Aa = 15.2±0.4MeV. The tangent of the 
scalar form factor at t = is known to be related to the scalar square radius. It is almost 
two times larger than the electric one, i.e. the (r 2 )fj ~ 1.6fm 2 while (r 2 )ff ~ 0.74fm 2 . The 
prediction of our model for the (r 2 )ff is 1.5fm 2 which is almost the same as obtained by Gasser 
and Leutwyler. It implies that the tail of the scalar density is of great importance. In Fig. 
3 we can find a long-stretched and strong tail in the sea contribution to the scalar density. 
This tail is due to the mesonic clouds arising from the Dirac sea polarization. Moreover, 
the sea contribution in the scalar density is large, compared with the other densities such 
as electromagnetic densities p3|-p5||. 

The other interesting quantities are presented in table 1. <Jq is the condensate of the 
singlet scalar quark operator in the nucleon:a"o = rh(N\uu + dd + ss\N) R s is defined by 
R s = (N\ss\N)/(N\uu + dd + ss\N). 
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IV. SUMMARY AND CONCLUSION 



We have discussed the scalar form factor with related quantities in the SU(3) NJL soliton 
model. The results we have obtained are in a good agreement with empirical data |3,21 



The reliable strangeness contents of the nucleon in the scalar channel is obtained by taking 
into account the m s corrections, since they suppress the excess of (iV|ss|iV) due to the 
leading order and rotational l/N c contributions. In contrast to Refs. ||0 suggesting no 
strangeness contribution, our model favors y = 0.27. The large value of the (r 2 )fj is caused 
by the pronounced long ranging tail which can be identified with the pion and kaon clouds. 
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TABLES 

TABLE I. The physical quantities related to the scalar form factor. The empirical data come 
from Ref.[3,17]. 



M 


370 MeV 


420 MeV 


450 MeV 


H/Xp 


fil s [IVie V J 





156.75 





148.49 





145.35 


rr .-r\\Ra\.T\ 

v-kN [ivievj 


A o no 


A A ^71 

44. 1 1 


/in ni 
40.01 


A n on 
40.80 


qo on 


38.69 


410 it o 


a [MeV] 


53.25 


49.25 


49.58 


46.24 


47.37 


44.35 




a g [MeV] 


22.77 


35.63 


20.87 


29.92 


19.92 


28.37 




y 


0.47 


0.20 


0.48 


0.27 


0.48 


0.29 


0.2 ±0.2 


Rs 


0.19 


0.09 


0.19 


0.12 


0.19 


0.13 




Act [MeV] 


32.29 


33.37 


18.36 


18.18 


14.23 


13.84 


15.2 ±0.4 


(r 2 ) S N 


1.94 


1.87 


1.56 


1.50 


1.40 


1.34 


1.6 
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Fig. 1 The scalar form factors of the nucleoli cr(t): The dotted curve corresponds to the 
constituent quark mass M = 370MeV, while the solid curve is for M = 420MeV. The 
dashed curve displays the case of M = 450MeV. The error bar denotes the empirical value 
from Gasser et al. [3]. 
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Fig. 2 The scalar form factors of the nucleoli cr{t): The solid curve displays the case 
m s = 148.49 MeV, while the dashed curve shows the case of m s = MeV. The error b 
denotes the empirical value from Gasser et al. [3]. 




Fig. 3 The scalar densities for the scalar form factor: The solid curve draws the valence 
contribution, whereas the dashed curve shows the sea contribution appeared from the 
mesonic clouds, i.e. Dirac sea polarization. 



